Electrical Detection of Spin Excitations 
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The objective of this paper is to draw attention to a possible new approach for measuring the 
spin excitation spectrum of a spin array by placing it in intimate contact with the channel of a 
spin-valve device in an anti-parallel (AP) configuration. We show, using realistic device parameters, 
that the second derivative of the current d^I /dV^ should exhibit peaks whose location corresponds 
to the energy of excitation and whose height is its strength. The effect is maximum for ideal half 
metallic contacts in AP configuration but we show that even with contacts that are less than ideal, 
the spin signals can be extracted using a de-embedding scheme. Compared to existing techniques 
our proposal has three attractive features: (a) high sensitivity down to single spins, (b) high spatial 
resolution limited only by the wavefunctions of the individual spins and (c) a relatively common 
device structure to which a lot of current research is devoted 
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PACS numbers: 72.25.Dc 

In recent years there has been increased interest in the 
dynamics of spin arrays, partly motivated by the possibil- 
ity of using such arrays for a wide range of applications, 
such as quantum computingiiSiii, data storage^, biolog- 
ical sensing^, DNA tagging and targeted drug deliverjii. 
Except for data storage, these applications require the 
ability to detect the spin spectra of only a few spins. 
This represents a major challenge^ and there has been 
a number of new approaches which have attempted 
to improve on the standard techniques of spin detec- 
tion which rely on detecting electronic or nuclear spin 
resonanceiSiiiiiSiiiiiii^. The current state of the art is 
magnetic resonance force microscopy in which single spin 
resonance detection^ has been recently reported. An al- 
ternative approach is to couple the resonance to the cur- 
rent flow in an adjoining conductor. Engel and Lossii 
proposed a method to detect the single spin decoherence 
inside a quantum dot by looking at the charge current. 
More recently, Xiao et. aM reported detection of sin- 
gle spin resonance by looking at the current that flows 
through a conventional Si MOSFET. 

However, a common problem of any type of resonance 
detection is the spatial resolution^fi. With the standard 
techniques, it is very difficult to achieve a resolution 
below Ifimm. Even with magnetic force resonance mi- 
croscopy, which was originally proposed to improve the 
spatial resolution, the state of the art is 25 nm. The 
lack of spatial resolution a well-known matter of concern 
in detecting very small magnetic particles using MRli^. 
For quantum computing algorithms too, one often needs 
to be able to readout a specific portion of the spin array 
with high spatial resolution. 

In this paper, we show that a lateral conductor having 
spin-polarized contacts in AP configuration (see Fig. 
allows one to detect a very small (even one) number of 
spins with a spatial resolution limited only by the ex- 
tent of the wave function of the individual spins. Ever 
since the earliest proposal for a "spin-transistor"iii, de- 
vices with ferromagnetic contacts have been the focus of 
much research effort. However, most of the work is de- 



voted towards making parallel contacts. We show that 
such parallel configurations are unable to detect the spin 
spectra, which should show up if we use the AP configura- 
tion instead. In this configuration, the second derivative 
of the current with respect to voltage (d^I/dV^) provides 
a direct measurement of the spin excitation spectrum. 
However, if the contacts do not have perfect injection 
efficiency, the spin signatures may get mixed up with 
contributions from other effects e.g. electron-phonon in- 
teraction. For this case, we propose a 'de-embedding' 
technique that can recover the spin information. 

FigQlshows the proposed device structure which could 
be a spin- valve device^ or a commercial FET with ferro- 
magnetic( ideally half-metallic) contacts. The spin array 
could, for example, be the paramagnetic traps at the Si- 
Si02 interface of an FET—, ferromagnetic atoms trapped 
inside a host molecule attached to the channeli^, or a fer- 
romagnetic gate with good exchange coupling with the 
channel electrons. 

Ideal Contact. - First we assume ideal half metallic con- 
tacts such that the source injects (and detects) only '-|-z'- 
polarized spins, while the drain injects (and detects) only 
'-z'-polarized spins. The current can flow only if elec- 
trons can flip their spins by interacting with a spin ar- 
ray.The point is that in this idealized structure any cur- 
rent that flows is a measure of the spin-flip scattering 
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FIG. 1: A schematic showing the top view of the proposed de- 
vice. The channel is connected to ferromagnetic/halfmetallic 
leads arranged in an anti-parallel configuration. A spin ar- 
ray, is placed in close proximity of the channel and affects the 
conduction by spin flip scattering through exchange coupling 
with channel electrons. 
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in the channel region and as such provides a measure of 
the excitations of the spin array. In our model, we as- 
sume the spin-scattering is caused by an exchange inter- 
action between a channel electron and individual spins 
of the array. The form of the exchange interaction is 
taken to be X]_r ~ ^j)^ ' where, r and Rj arc 
the spatial coordinates and a and Sj are the spin opera- 
tors for the channel electron and j-th spin in the array. 
J{f^Rj) is the interaction constant between the channel 
electron and the j-thc spin in the array. The current flow 
in the device is proportional to the difference between the 
rates at which the channel electrons are scattered from 
'-z'^'+z'(dcnoted by .Fi) and '+z'^'-z'(dcnoted by F2): 
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Is = e— {Fi ~ F2) 



(1) 



where Fi and f 2 are given by 

Fi= \<m\,(3\J{f^Rj)3-Sj\ni,a>\'^ 

a,l3,'m\,ni,Rj 

5{Ej3 -f — Ea~ Enl)Pafnl{i — /jrit) 

(2a) 



a,p..m] ,ni.,R 



< n l,a\J{f ~ Rj)d ■ Sj\m 13 > I'' 

S{Ea + Enl — Efj ~ Em-\)Pflfm\{^ — fnl) 



Here, a or /? is an eigen state of the spin array with 
occupation probability Pa such that = 1; rn t 

and n J, describe the '+z' and z' eigcnstates of the 
electron system with occupation functions f„i^ and /„| 
respectively. Substituting a.Sj = OzSz.j + 1/2(ct+5_j + 
cr_S'+j ) and replacing sums over m, n with integrals over 
the density of states D{E), we obtain 



Fi = j dE Jo^DiiE) J d{hLj)D^{E ~ huj) 
fl{E){l-f^{E-hw)}A^{u:) 

F2 = j dE Ja^Di{E) J d{haj)D^ {E - hu) 

f^{E-nuj){l-f^{E)}A+{^uj) 
where we have defined: 



(3a) 



(3b) 



A_(a;) =^F(i?,)^F„l </3|5_,,|a> pJ(£;^-F,-M 

Rj a,l3 

(4a) 



Here, 1/4| / V:,(r) J(r-^,)^„(f)drp = J^FiR,), Jq be- 
ing a a constant and F{Rj) representing the functional 
variation of interaction with Rj. Strictly speaking, F{Rj) 
can depend on m,n; but we have ignored this in the fol- 
lowing treatment. Clearly A_(tt') and A^{~uj) describe 
the spin excitation spectra of the spin array. Fig. I3a),(b) 
show the current vs. voltage calculated from Ea. Hll3|l for 
an array of Ni spins having excitation spectra Az^(±uj) 
as indicated in the figure. Note that the Azp{±uj) in case 
1 corresponds to an unpolarized (non-interacting) array, 
while those in cases 2 and 3 correspond to arrays polar- 
ized in '-z' and 'z' directions respectively. Unlike case 1, 
which leads to a linear I-V, cases 2 and 3 show rectifi- 
cation because a specific polarization allows current to 
flow by spin flip scattering only for injection from one 
contact but not the other. The curves in Fig.(|2Ia)) as- 
sume Ea = Ep, whereas for the curves in Fig.(|2fb)). 
Ef3 — Ea = hw = SkT, resulting in a simple shift in the 
I-V characteristics along the voltage axis by 3kT. 

For the plots in Fig. 2, we have assumed that A±{'^uj) 
remains unchanged despite current flow. This will be true 
if the relaxation processes (say denoted by 7/ at w) are 
strong enough to keep A±{^uj) close to their equilibrium 
values Aj.{^uj) despite the spin-flip scattering process 
induced by the current. By writing the rate equations for 



A_(w) and A+^—Lo) and noting that A^{uj) - 
a conserved quantity one can show thatiiS 



'A, 



-Lo) is 



(2b) A^iu;) - A+i-Lo) 



(7+ - 7-) 
7+ + 7- + 7/ 
7/ 

7+ + 7- + 7/ 



{A^iuj) + A+i-Lu)} 

{A"_ico) - Ali-u;)} . 

(5) 



Here, 7_ and 74. arc found respectively from Eas. l|3a|) 
and Ij3bp by excluding (i)thc integral over d{huj) and (ii) 
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1 : A_(a) = A^{-a) = 0.5N,S(tico) 
2 : A^(-co) = N,d(ho)); A_(co) = 
3 : A _(co) = N,d(nco); A^ (-co) = 
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4 : A_{m) = A^{-m) = 0.5N,S(ti(O-3kT) 
5 : = e-'""''^A_(0)) = N,S(tico-3kT) 
6 : A_(oi) = e-'"'"*^4^(-w) = N,S{tico-3kT) 



FIG. 2: I-V characteristics for different spin excitations ac- 
cording to Eg. (1113 bH for ideal contact. The simulation is done 
using the following parameters: D = m* LW/{2-Kh^), where 
L = lOOnm and W — 300 nm, while m* = 0.2mo, mo being 
A+{-uj) ^ F{Rj) Y Pf^l < a\S+^j\p > \^5{Ep-Ea-hi the free electron mass. Ni = NiLW where ffi « lO^/cm^, 
ji- a p T = 300K and Jo = 0.02/ieV which corresponds to an ex- 

(4b) change frequenc}«i of 30 MHz. 
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A-{uj) (for 7_) and A+(-uj) (for 7+). From Eq.Q, 
A ± (=Fw) will remain unaffected by the current if 77 >> 
7+ + 7- . For the parameters used in Fig. ^ (see caption) , 
this requires 7/ to exceed {h/10 ns) (T/SOOK) which is 
quite likely in practical systems. Note that, in this case, 
the currents shown in Fig.(|21) are steady state currents. 
But if the condition 7/ >> 7+ + 7_ is not fulfilled, the 
currents should be viewed as the transient currents that 
flow initially at t=0. Indeed, this effect can then be used 
to write information into the spin array^SiSi. 

The integrands of Fi and F2 in Eqs.(|3J) cancel each 
other at all energies if 



A+{-oj) \l-f^ 



h 



(6) 



where = f^{E — huj). This means that the current for 
zero bias is zero (see Eq.(^) if the spins are in equilibrium 
with the same temperature so that yl_([j)/A+(— = 
e'!cp{huj / kT) . But if external forces (say, for example, a 
resonant r.f. field), drive the spins out of equilibrium, 
then the zero-bias current need not be zero and under 
the right conditions it may be possible to extract energy 
from the spins and deliver it to an electrical load. For 
cases 2,3,4 and 6 the spin array is not in equilibrium with 
the thermal bath, giving rise to a non-zero current at zero 
bias. 

Non-Ideal Contact.- Real systems are likely to have 
limited spin injection efficiency, resulting in a ballis- 
tic current in addition to the scattering current. This 
calls for a generalization of the expression for current 
(Eq.(^)). If we define the contact couplings for '±z' spins 
by Ti^/n-^/i^, it can be shown using the Non-Equilibrium 
Green's Function (NEGF)method that the terminal cur- 
rent is given by 



dE 



(7) 



where /| and fi have to be determined sclf-consistcntly 
from a combination of rate equations describing in and 
outflows and current conservation within the device^. 
From Eq.(|Zj), it can be shown that 



(8) 



where 

Ib = 

and 



-r / dE 
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Ul - Ir) 



dE 
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Eli + Ffii Tli + Tri 



(9) 
(10) 



In Eq.^,is{E) is defined as e(27r/ri) |Fi(£') - F2(£')|, 

where Fi{E), and F2{E) arc the integrands of the vari- 
able 'E' (terms inside the '[ ]') in Ea. l|lla|) and (|3H|) re- 
spectively. The first term in Eq. (|SJ) is the ballistic current 
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FIG. 3: I-V calculated for Non-Ideal contact from Eq.©. 
The injection efficiencies of the source and drain were as- 
sumed to be 60% and 40% respectively. Two spin excita- 
tions of hui — 2.5 and 5 meV were assumed. (a)Total cur- 
rent. No features can be detected due to large contribution 
from ballistic current, (b) The second derivative of the total 
current. The peaks retain the information of spin scatter- 
ing, (c) dotted line — anti-symmetric part of the total cur- 
rent: Jasym = I (V) — I{—V); solid line — second derivative of 
the anti-symmetric part, (d) dotted line — symmetric part of 
the total current: /sym = I{V) -\- I{ — V); solid line — second 
derivative of the symmetric part. The contact couplings are 
roughly around 0.1 eV. All other transport parameters are 
the same as were indicated in Fig|5| The temperature was 
set to T=3K. This is necessary to distinguish among the mul- 
tiple peaks which will be otherwise washed away because of 
temperature broadening (see Eg. 1141 '). 



(Ib) which is clearly a sum of the currents that flow from 
left to right in the 'up-up' and the 'down-down' channels. 
On the other hand, the second term l^ results from the 
scattering current, /^(E') flowing between the 'up-down' 
channels and gets distributed according to the conduc- 
tance ratios of left and right. 

d^I/dV^ as a signature of spin excitation.-Vnder small 
voltage bias (this is the case for all of our treatments), 
Ib is linear in applied bias,!^, whereas 7^., having terms 
containing product of two fermi functions, is quadratic in 
V. Therefore, even though all features coming from spin 
scattering in the current itself is essentially hidden by 
a more dominant ballistic part (see Fig.Qa)), the second 
derivative very effectively retains the signature of the spin 
excitation (see Fig. |2Ib)). For simplicity, we assume 
the rate of spin flip scattering, which is proportional to 
j'^Azp{±uj), is small compared to Ti^^r.-^/^, so that w 
ai/l + bifR and fi ~ a2fL + &2/i?,, where 



ai 



02 



-; bi 

Flt + r^T + Fflj- 

Tii . , ^Ri 

02 



(11) 



Tlx + ^Ri ^Li + ^Ri 

We also assume that the density of states is constant in 
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the energy range of interest, which should normally be 
true for low bias. Taking the second derivative of the 
terminal current (Eq.®), noting that the ballistic part 
is purely linear and recognizing that the voltage enters 
Fi and F2 only through the fermi functions, we obtain 

0=e^ (y)^'^o/ d(M [A+{-u:)-A^{u:)] 
{aib2S{hiu + eV) + a2bi6{eV - hu)} 

(12) 

at T = 0. Note that (From Eq.©) 

[A+i-u;) - A„(^)] = J2 ^(^^■) E I < ^l^-l" > I' 
5{Ep-E^-huj)[Pfi~P^) 

(13) 

is a measure of the spin excitation spectrum and is re- 
flected in the amplitudes of the two delta functions at 
eV = ihio, only one of which is non-zero for ideal con- 
tacts. This is the main result of this paper. When T ^ Q 
the S is broadened into a convolution of two sech^ func- 
tions: 

/.Esoc„'(^).oci^ (51^), 

(14) 

which describes the numerically computed line shapes in 
FigEl quite well. 

Symmetrization- In the above derivation, we have 
dropped d'^IsldV^ arguing that Is vs. V represents 
the current in the absence of spin-flip scattering and 
is exactly linear (since D is assumed to be constant). 
But, if scattering processes that do not flip spin (like in- 
elastic phonon scattering) are present, d'^Is/dV'^ is not 
zero. However, it is antisymmetric in V (See RcfsiSSiS) 
and therefore, can be eliminated if we take the sym- 
metric component of the second derivative by writing 
d^hy^/dV^ = d^jdV^ {liy) /(-^)}, so that the con- 
tribution from d^lB/dV^ cancels out. 

Discussion.- Note that the current calculated in Fig. 
2 arises from Nj = 30 spins in the array interacting with 
the channel. This current is many orders of magnitude 
larger than the measurement limit, even though we have 
used a very conservative estimate for the interaction pa- 
rameter Jq. This means that even single spins should be 
detectable in realistic structures. 

We believe that a significant strength of the proposed 
techniques stems from its potential spatial resolution. 
Since the interaction constant J is determined by the 
wavefunction overlap of the spins with the channel elec- 
trons, we expect the spatial resolution of our proposed 
method to be limited only by the spread of the wave- 
functions of the individual spins, which is the ultimate 
limit on any measurement technique. For shallow donors 
the spin wavefunction extend over distances ~ (Bohr 



radius) which is lOOA" in GaAs. Deep level traps have 
even more localized wavefunctions. 

One requirement for our proposed device is that the 
channel length should be much less than the spin- 
relaxation length due to other causes like spin-orbit cou- 
pling. However, this poses no serious restrictions. For 
example, in doped GaAs, spin relaxation lengths over 1 
/im has been reported24. Other channel materials of in- 
terest like Si, carbon nanotubes etc. are known to have 
much weaker spin-orbit coupling. 

Ideal half metallic contacts in the AP configuration 
leads to r^t = Tm ^ and F^j — Vn^ — 0. As seen 
from Eq.||Sl), this situation maximizes the effect of /, on 
the terminal current (/^ = J dEIs{E); Ib = 0). If, on 
the contrary, Tli/Tli = Tji-^/Tj^i, Is makes no con- 
tribution to the terminal current (/, = 0). This cor- 
responds to symmetric unpolarizcd contacts, which com- 
pletely destroys the effects we are seeking to exploit. This 
is in qualitative agreement with a previous resultS^ that 
showed, in connection with the 'zero bias anomaly' in 
tunneling magneto resistance (TMR) devices, that spin- 
flip scattering is more prominent when the spin-polarized 
contacts are anti-parallel (AP) rather than parallel(P). 

We emphasize the fact that the use of the ferromag- 
netic contacts make it possible to detect the spin spec- 
tra without having to put an external microwave signal. 
However, the proposed device could be combined with a 
microwave signal to measure the spin decoherencc time 
T2. Near the resonance condition, the total deviation 
of the average spin from the direction of static mag- 
netic field depends on the frequency of the rf field, Wrf^. 
As a result the current as a function of cuy-f will show a 
dip whose width is proportional to the transverse deco- 
herencc time T2. Note from the example in Fig. 2, the 
changes in the current levels is in the order of nA, which 
is easily detectable^. 

Our basic result for the current (see Eq.®) can be 
understood by combining a golden rule treatment of 
the spin-flip current (see Eq.lQ)) with a heuristic rule 
for distributing it between the contacts (see Ea. (|10() '). 
Note, however, that this result is obtained rigorously 
from a NEGF treatment. We have also obtained simi- 
lar results directly from a Landauer-type scattering for- 
mulation with entangled 'channel electron'- 'spin array' 
states^. 

In summary, we have proposed a technique which 
should allow one to detect spin excitation spectra of a 
spin array. Compared to existing techniques our proposal 
has three attractive features: (a) high sensitivity down 
to single spins, (b) high spatial resolution limited only by 
the wavefunctions of the individual spins and (c) a rela- 
tively common device structure to which a lot of current 
research is devoted^!. Possible applications include sin- 
gle spin detection, relaxation time measurements, biolog- 
ical sensing of miniature magnetic particles and read out 
devices for quantum computing, high density memory^, 
spin-wave interconnect and logic applications^^. 
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See Eq.(4) in—. The rate equations and resulting self- 
consistency were discussed for the special case where '±z' 
states are degenerate. The subscripts 'u' and 'd' in corre- 
spond to our I and \ respectively. 

a longer version showing these comparisons will be pub- 
lished elsewhere. 



